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Abstract 



We demonstrate that measurement of tagged structure functions of the 



\Q . deuteron in (e, e'N) semi-inclusive reactions can discriminate between dif- 

ferent hypotheses on the origin of the nuclear EMC effect. By choosing ex- 
treme backward kinematics for the spectator nucleon to minimize effects from 
the deuteron wave function and final state interactions, one can isolate the 
modifications in the structure of the bound nucleon within the impulse ap- 
proximation. The same reaction can be used to extract the large- x neutron 
to proton structure function ratio. 

X 

I. INTRODUCTION 

More than a decade after the discovery of the nuclear EMC effect |J and many fine mea- 
surements of the ratios of structure functions of nuclei and the deuteron, no consensus 
has been reached on the origin of the effect. The x dependence of the effect, while non-trivial, 
is rather smooth and has the same basic shape for all nuclei, making it is easy to fit in a 
wide range of models with very different underlying assumptions. The only extra constraint 
available so far comes for measurements of the A-dependence of the sea distribution, which 
restricts some of the models, but is still not sufficient to allow one to unambiguously identify 
the origin of the EMC effect. 

In order to move beyond this rather unsatisfactory situation, new experiments involving 
more kinematical variables accessible to accurate measurements are necessary. The aim of 
this study is to demonstrate that use of semi-inclusive processes off the deuteron, 

e + D ^e + N + X, (1) 
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where a nucleoli is detected in the target deuteron fragmentation region, may help to dis- 
criminate between some classes of models. In particular, one may be able to distinguish 
between models in which the effect arises entirely from hadronic degrees of freedom — nu- 
cleons and pions (which in the traditional nuclear physics picture are responsible for the 
binding of the nucleus), and models in which the effect is attributed to the explicit defor- 
mation of the wave function of the bound nucleon itself. By selecting recoil nucleons with 
small transverse momentum in the backward region, effects due to final state interactions 
(FSI) can be minimized, thus allowing one to probe, within the impulse approximation, the 
deformation of the structure of the bound nucleon. 

It is worth emphasizing that measurement of the reaction ([!]) in the kinematics of interest 
will become feasible in the near future, for example at Jefferson Lab ||. Recent development 
of silicon detectors also make it possible in the jet target experiments to measure recoil 
nucleons with a low momentum threshold of about 100-150 MeV/c. Plans to implement 
such a technique are under discussion at the HERMES detector at HERA J7J. 

Aside from providing insight into the origin of the nuclear EMC effect per se, the mea- 
surements of tagged events may also be useful in connection with the problem of extracting 
the neutron structure function from deuteron data. By selecting only the slowest recoil pro- 
tons in the target fragmentation region, one should be able to isolate the situation whereby 
the virtual photon scatters from a nearly on-shell neutron in the deuteron. In this way one 
may hope to extract the F 2n structure function while minimizing uncertainties arising from 
modeling the nuclear effects in the deuteron. 

This paper is organized as follows. In Section II we outline the basic formalism for semi- 
inclusive deep-inelastic scattering off the deuteron, and discuss the conditions under which 
the impulse approximation may be valid. Section III is devoted to a survey of several models 
whose predictions for the ratios of tagged structure functions are compared. Possibilities of 
extracting the neutron structure function from semi-inclusive experiments are discussed in 
Section IV, and finally some conclusions are drawn in Section V. 



II. BASIC FORMALISM 

It was observed a long time ago that the nuclear EMC effect for the deviation from unity 
of the ratio 

is approximately proportional to the nuclear density. This is natural for a dilute system, 
and indicates that most of the EMC effect is due to two-nucleon interactions. Based on 
this observation one should expect that the EMC effect in the deuteron is much smaller 
than that in heavy nuclei. However, by virtue of the uncertainty principle, one may try to 
enhance the effect by isolating the configuration where the two nucleons in the deuteron are 
close together. For example, it is easy to check that the main contributions to the deuteron 
wave function for nucleon momenta p > 300 MeV/c come from distances < 1.2 fm. 

To make our discussion quantitative, let us begin by writing the electromagnetic tensor 
Wp U of the deuteron in terms of the matrix element of the electromagnetic current: 
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= E I I X,N)(X 7 N | Jtf(0) | D). (3) 

The tensor can be expanded in terms of four possible Lorentz structures, with the 
coefficients given by the invariant structure functions F®, F®, F® L and F® T (T = transverse, 
L = longitudinal). These functions can in general depend on four variables, constructed from 
the four- momenta of the target deuteron (P), virtual photon (q) and spectator nucleon (p s ) 
(or equivalently the momentum of the struck nucleon p, where p = — p s in the deuteron rest 
frame). In terms of the invariant structure functions the differential cross section for the 
semi-inclusive reaction (e, e'N) can then be written: 



da 2a 2 em ( i _ y _xWM 2 " 



dxdQ 2 d 3 p s /E s xQ 4 \ Q 2 



FE + ( ^ + tan 2 (£/2) 1 —F* + ( ^ + tan 2 (0/2)J cos0 Ff L + cos(20) F? T 



(4) 



where the structure functions are related to the components of the electromagnetic tensor 
W£ v by: 

F£ = + cos<5) 2 ■ WE~, (5a) 
F$ = M{Wf) X + W y D v ), (5b) 
F£ L = 2z/(l + cos 5) ■ Wo x \ (5c) 
F° T = ~ sin 2 5 ■ (W^ - W y D y ) . (5d) 

The kinematic variables in Eqs.(^j) and (^) are x = Q 2 /2Mu, where v is the energy of virtual 
photon in the target rest frame and M the nucleon mass, Q 2 = 4E(E — v) sin 2 (6 l /2) is the 
squared four-momentum transfer to the target, E is the beam energy and y = vjE. The 



angle (p is the azimuthal angle for the spectator nucleon, and E s = y M 2 + p 2 the spectator 

nucleon energy, and sin 2 5 = Q 2 /q 2 in Eqs.(|5]). We define the photon three-momentum q to 
be in the +z direction. 

A convenient choice of the four independent variables for the structure functions is the 
two inclusive deep-inelastic scattering variables, x and Q 2 , and the transverse momentum, 
Py, and light-cone momentum fraction, a s , of the spectator nucleon: 

E.s-P s z , R , 
«. = (6) 

where p s z is the longitudinal momentum of the detected nucleon (to simplify the expressions 
we have neglected here the deuteron binding energy, Mo « 2M). In the Bjorken limit the 
variable a s then satisfies the condition ||: 

a s <2-x. (7) 

Having defined the relevant cross sections, we must now establish the kinematical range 
in which the nuclear modifications of the bound nucleon structure function will be accessible 
experimentally. From Eqs.(H) and (Bh the restrictions: 
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g 2 /q 2 ~ 4m 2 x/q 2 < 1, 

Q 2 < E(E-u), 



(8a) 
(8b) 



enhance the contribution of the longitudinal structure function Ft, which is expressed 
through the "good" component of the electromagnetic current, Eq.(^). Another simpli- 
fication can be achieved by considering the situation where the detected nucleon in the 
deuteron is in the spectator kinematics, namely: 

a s >l-x. (9) 

In this kinematical region the contribution of the direct process where a nucleon is produced 
at the 7*iV interaction vertex is negligible ||. 



A. Factorization and the Impulse Approximation 



In formulating deep-inelastic scattering from the deuteron the simplest approach adopted 
has been the impulse approximation for the nuclear system. We will consider two formu- 
lations of the impulse approximation: one based on the covariant Feynman, or instant- 
form, approach, where one nucleon is on-mass-shell and one off- mass-shell, and the non- 
covariant/Hamiltonian light-cone approach, in which both nucleons are on- mass-shell (but 
off the light-cone energy shell). In the next subsection we shall consider corrections to the 
impulse approximation, in the form of final state interactions, but for now let us review 
briefly the basic impulse approximation assumptions and results. 



1. Covariant Instant- Form Approach 

In order to construct covariant amplitudes in the instant-form of quantization, summa- 
tion over all possible time-orderings of intermediate states is essential. Incorporation of 
negative energy configurations into the total Lorentz-invariant amplitude is done on the 
basis of introducing intermediate state particles which are off their mass shells. In the 
impulse approximation for the scattering of a virtual photon, 7*, from the deuteron, the 
invariant amplitude for the complete process factorizes into a product of amplitudes for 7*- 
off-mass-shell nucleon scattering, and for forward nucleon-deuteron scattering. In this case 
the hadronic tensor of the deuteron can be written || : 

W£ u (P,p,q) = Tr [S ND (P,p) W N u (p,q)} , (10) 

where the nucleon-deuteron scattering amplitude Snd (for spin-averaged processes) contains 
scalar and vector components: 

S ND = S + 7^. (11) 

The operator W N U in Eq. ( |ToD is the truncated hadronic tensor for the off-shell nucleon, which 
describes the 7* A" interaction. Because W N V is a matrix in both Lorentz and Dirac spaces, 
its structure is necessarily more general than that for a free nucleon. In the Bjorken limit, it 
was shown in Ref. IJ that out of a possible 13, only three independent structures contribute: 
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Wj?(p,q) = -(g^ + ^ (Wo(p,q)+^W 1 (p,q)+^W 2 (p,q)) + O (J^J . (12) 

Note that the full expression for W^ u contains both positive and negative energy pieces. 
The on-shell nucleon tensor is obtained from W^ u by projecting out the positive energy 
components: 

W^{p,q) = ]Tr \(tf + M)W^(p,q)} (13a) 
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9^ + ^ j (MW + M 2 Wi + p ■ qW 2 ) , (13b) 



where the functions W0-2 here are evaluated at their on-shell points. 

In terms of the off-shell truncated functions W0...2 the deuteron tensor can then be writ- 
ten: 

Wj?(P,P, q) = -4 (g^ + ^ (So W Q + S ■ p W x + S ■ q W 2 ) , (14) 

so that in general the right-hand-side does not factorize into a single term with separate 
nuclear and nucleon components. As explained in Ref. , factorization can be recovered by 
projecting from «So.i only the positive energy (on-shell) components, in which case: 

So oc S ■ p oc S ■ q. (15) 

This proportionality can also be obtained by taking the non-relativistic limit for the ND 
amplitudes to order p 2 /M 2 |TIJ, since all negative energy contributions enter at higher orders 
ll]. 12]. For the relativistic ND interaction, the negative energy components enter through 



the P-state wave functions of the deuteron. However, in most realistic calculations of these, 



in the context of relativistic meson-exchange models of the deuteron [13], the contribution 
of the P-state wave functions is only a fraction of a percent, and so long as one avoids 
regions of extreme kinematics which are sensitive to the very large momentum components 
of the deuteron wave function, dropping the negative energy components is a reasonable 
approximation. Indeed, the explicit calculations of the inclusive structure function of the 
deuteron with relativistic wave functions shows that the factorization-breaking corrections 
amount to < 1-2% for all values of x ^5 0.9 [11,13]. 



Finally, given these approximations one can write the deuteron hadronic tensor as: 

W^(P,p,q) « S IF (P,p) eff (p,q), (16) 
where S IF is the nucleon spectral function within the instant-form impulse approximation 
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S IF (a s ,p T ) = (2 - *s) 2{ ^_ D Es) \M*s,Pt)\ 2 , (17) 



normalized such that 
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/d 2 p r — S IF (a a ,pT) = l 
J a., 



(18) 



to ensure baryon number conservation ^5||16j. Note that the approximation for the spectral 
function in Eq.fllTD is valid only for non-relativistic momenta, and for large momenta the 
full (non-convolution) expression for the hadronic tensor in Ref. (T^j should be used. The 
effective nucleon hadronic tensor Wj? e ^ is defined as: 



= - y9nu + 



(MWo + p 2 W x + p-qW 2 



P eff 



X 



2-«, 



,P Z ,Q Z + 



(19) 



where the effective nucleon structure function F^f is now a function of the momentum 



fraction xj (2 — a s ) and the virtuality p of the bound nucleon, as well as Q . The kinematics 
of the spectator process gives rise to the relation: 



V 



2p\ + (2 - a s )M 2 



+ -(2-a s )M\ 



D J 



(20) 



using which one can equivalently express F^l/ as a function of the transverse momentum px 
of the interacting nucleon, rather than p 2 . 



2. Light-Cone Approach 

In the light-cone formalism one can formally avoid the problems associated with negative 
energy solutions in Eq.(|l"2"D, however, to obtain Eq.flTo]) one must on the other hand consider 
contributions arising from instantaneous interactions [17||. To take these effects into account 



one has to use gauge invariance to express the contribution of the "bad" current components 
of the electromagnetic tensor through the "good" components: 

Ji = ~^ J t (21) 

and include the contribution of the instantaneous exchanges using the prescription of Brod- 
sky and Lepage 0] . In the approximation when other than two-nucleon degrees of freedom 
in the deuteron wave function can be neglected, one can unambiguously relate the light-cone 
deuteron wave functions to those calculated in the rest frame in non-relativistic instant-form 
calculations [|1£],^(|. The final result for the spectral function in the light-cone approach can 
be written similarly to Eq . (|T6l) , only now S is replaced by the light-cone density matrix: 



where 



S LC (a s , PT ) = ^-^\Mk)\\ (22) 
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^±it _ AP (23) 

\ a s {2 - a s ) 



is the relative momentum of the two nucleons on the light-cone, and eff in Eq.© 



is now the bound nucleon electromagnetic tensor defined on the light-cone. Note that the 
main operational difference between the instant-form (fl7|) and light-cone ( p2|) impulse ap- 
proximations is the different relation between the deuteron wave function and the scattering 
amplitude. Numerical studies have demonstrated that in the kinematical region of interest, 
|p a | ^5 0.5 GeV/c, the difference between the results of the two approximations (for the same 



deuteron wave function) is quite small — see Sec. |III D| and Ref. || . 

The structure functions for the scattering from the off- "+" -shell bound nucleons may 
depend on the variables of this nucleon similarly to the case of the Bethe-Salpeter, or co- 
variant Feynman, approach. In another language this dependence can be interpreted as the 
presence of non-nucleonic degrees of freedom in the deuteron. With this is mind, we shall 
use Eq. fll6|) as the basis for the results discussed in the following sections. Before focusing 
on specific model calculations of the semi-inclusive deuteron structure functions, however, 
let us first turn our attention to the validity of the impulse approximation, and the problem 
of final state interactions in particular. 

B. Final State Interactions 

In the kinematical region defined by Eq.(^) the contribution of direct processes, where a 
nucleon is produced in the 7*iV interaction, is negligible ||. Therefore within the framework 
of the distorted wave impulse approximation (DWIA) [] the total deuteron tensor Wj^ can 
be expressed through the nucleon electromagnetic currents as: 

W£ u (x jasjPTj Q 2 ) w \j2(D\pn)(XN | O ia + O FS i | XiV)| 2 Wj^(x, a s ,pT, Q 2 ) 

= S DWIA (a s ,p T ) eff (x,a s ,p T ,Q 2 ), (24) 

where Oia is the impulse approximation operator, while Op si describes the soft final state 
interactions between the final hadronic products and the spectator nucleon. The function 
S DW (a s ,pT) now represents the spectral function distorted by FSI effects. 

Analysis |2T] of the recent high energy deep-inelastic scattering data on slow neutron 
production |22| is rather indicative that even in heavy nuclei final state interactions are 



small, so that the average number of hadrons which reinteract in the target does not exceed 
unity. This indicates that the system which is produced in the 7*iV interaction is quite 
coherent and interacts at high energies with a relatively small effective cross section, 

(Jeff < &NN- (25) 



1 Note that the DWIA approach works best at extreme backward kinematics, Eq.(p6|) below, where 
final state interactions have a small contribution, Eq.(|27|). 
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Such a situation allows one to factorize the 7*iV interaction and use the calculation of FSI 
in eD — > e p n processes as a conservative upper limit. 

The final simplification which we can gain is to consider the extreme backward kinematics 
where, in addition to Eq.(||), we also require: 

p T « 0. (26) 

In this case it can be shown ||2~3lJ2"3II that: 

a eff (Q 2 ,x) \tp D (a s ,<p T >)tp D (a s ,0)\ 



S DWIA (a s ,p T ^0)~S(a s ,p T ^0) 



1 - 



57T < r > 

pn 



S(a s ,p T ~ 0)/JE s E s (< p\ > 



(27) 

where < r 2 n > is the average separation of the nucleons within the deuteron, E s is the 



spectator nucleon energy, and E s (< p\ >) = y M 2 + p s z 2 + < p\ > is the energy evaluated at 

the average transverse momentum < p\ > 1//2 ~ 200-300 MeV/c transferred for the hadronic 
soft interactions with effective cross section <J e ff- The steep momentum dependence of the 
deuteron wave function, \ip£,{a S) < px >)| <C \ipo(o: s ,pT ~ 0)|, ensures that FSI effects 
are suppressed in the extreme backward kinematics, in which case the original impulse 
approximation expression for S(o> s ,Pt) can be used rather than S DW {ql s -iPt)- Finally, 
expressing the electromagnetic tensor of the nucleon, W^ v , through the effective nucleon 
structure function F^jJ and F^pf, we can then write for the deuteron tensor: 



V-Q _\ ( _ , P-q_ \ 1 rpeff ( x 



+ ^ + V*J ^ fe PT ' Q ))■ (28) 

Defining F^ttltt t° the semi-inclusive structure functions in Eqs.fj^) with the spectral 
function factored, 

Fl,t,tl,tt — S(a s ,pT) Fl,t,tl,tt > (29) 

one can express the semi-inclusive nucleon functions in terms of the effective structure func- 
tions of the nucleon as: 

i^(z,a s ,p T ,Q 2 ) = -sin 2 5-^ F*g (j^Pt, Q 2 ) 

+ (l + cos5) 2 (a s + ^c^ ^(^^,Pt,Q 2 ), (30a) 

F?(x,a s ,p T ,Q 2 ) = 2F&' (-^—,p T ,Q 2 ) + (-^—, P t,Q 2 ) , (30b) 

\2 — a s J M 2 v \2 — a s J 

F? L (x,a s ,p T ,Q 2 ) = 2(1 + cos5)g L + ^a q ) (^—^Q 2 ) , (30c) 

F^a s ,p T ,Q 2 ) = ^%\n S J (^T'PTiQ 2 ) , (30d) 



S 



where 

«9 = M » ( 31a ) 

_ p-g l + cos<5 M 2 +p2 

j/ = — — = q a s + a a — — . (31b) 

M IHI 2 q 2a s M v ; 

Note that within impulse approximation the Callan- Gross relation between the F^jJ and 
structure functions is preserved. Therefore the experimental verification of this relation 
could serve as another way to identify FSI effects. 

Equations (|30"D and ([H]) show that at fixed x and Q 2 — > oo, when a q — > the longitudi- 
nal structure function, F^, does not depend explicitly on the transverse momentum of the 
nucleon. It contains pr dependence only in the argument of the bound nucleon structure 
functions, which arises from the possible nuclear modifications of the nucleon's parton dis- 
tributions. The above argument, and the fact that FSI of hadronic products with p x ^ 1 
GeV/c practically conserve a s , allows one to conclude that the FSI effect on F^ is minimal. 

On the other hand the functions F T , F TL and F TT do explicitly depend on the specta- 
tor transverse momentum and therefore the hadronic reinteractions in the final state may 
strongly affect these structure functions. Such a situation suggests that a separate study of 
the complete set of the structure functions will allow one to investigate the effects of final 
state interactions in the deep-inelastic (e, e'N) reactions. Note that for the production of 
spectators with px > 0, FSI effects are not likely to depend strongly on x for x > 0.1. This 
is because at x > 0.1 the essential longitudinal distances in deep- inelastic scattering are 
small. Therefore for these kinematics one expects that a e ff(Q 2 , x) m a e ff(Q 2 ). The effect of 
FSI will be investigated in the dedicated deep-inelastic scattering experiments at HERMES 
which will measure the A-dependence of the forward produced hadrons. 

These observations enable us to conclude that in the kinematic region defined by Eqs.@, 
(S) and (PUD, where the contribution of Fjf is enhanced and FSI effects are small, the 
differential cross section (Ht) can be written: 



do" 3 -*** 2a* ,S(a s ,p T ) 



(i-y)- 



dxdW 2 d(\oga s )d 2 p T Q 4 v ^ 2 - a 



2q 2 M 



2n ' ,:l/ l-y)S(a s ,p^F^ (-^—,p T ,Q 2 ) (32) 

\2 — a s J 



where we have made the transformation dQ 2 /x — ► dW 2 , with W 2 = —Q 2 + IMv + M 2 . 

Based on the expectation that FSI effects should not strongly depend on x, from Eq.(p2|) 
it may be advantageous to consider the ratio of cross sections relative to a given x, in the 
range 0.1-0.2, where the observed EMC effect in inclusive scattering is small [|H| , f2"5| l : 

da(x 1 ,a S} p T} Q 2 ) I da(x 2 ,a s ,p T ,Q 2 ) 



G(a 8 ,p T ,x 1 ,X2,Q x 



dxdW 2 d(log a s )d 2 pT / dxdW 2 d(log a s )d 2 pT 
F^ f (x 1 /(2-a s ),p T ,Q 2 ) 



(33) 



F^ f (x 2 /(2-a s ),p T ,Q 2 Y 

In our analysis we will consider only production of backward nucleons, a s > 1, to suppress 
contributions from the direct processes where a nucleon is produced in the 7*iV interaction 
vertex. The more liberal condition, Eq.(|9]), is in reality sufficient |§. 
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Note also that for heavier nuclei the FSI becomes much more important in the limit of 



large x. As was demonstrated in Ref. ||17|| , in this limit rescattering of hadrons produced 
in the elementary deep-inelastic scattering off the short-range correlation is dynamically 
enhanced, since the average value of the Bjorken- variable for this mechanism is m x, as 
opposed to x/(2 — a s ) for the spectator mechanism. In this sense the deuteron target 
provides the best way of looking for the EMC effect for bound nucleons. 

III. MODELS 

In this section we briefly summarize several models of the EMC effect which we use in 
our analysis and present their predictions for the tagged structure functions. The differences 
between the models stem from dynamical assumptions about the deformation of the bound 
nucleon wave functions, and from the fraction of the EMC effect attributed to non-baryonic 
(mesonic) degrees of freedom in nuclei — from the dominant part in some versions of the 
binding model, to the models where non-baryonic degrees of freedom play no role, as in the 
color screening or QCD radiation models. (Other models which have been used in studies 
of semi-inclusive DIS include the six-quark cluster models discussed in Refs. p6|p7| .) 



A. Binding Models 

One of the simplest of the early ideas proposed to explain the nuclear EMC effect was the 
nuclear binding model, in which the main features of the EMC effect could be understood 
in terms of conventional nuclear degrees of freedom — nucleons and pions — responsible for 



the binding in nuclei [p!7| , p8| -|35|] . Within the formalism of Sec. [II A 1| , the inclusive nuclear 
structure function in the EMC ratio, Eq.(Q), is expressed through a convolution of the 
nuclear spectral function and the structure function of the bound (off-shell) nucleon (c.f. 



Eq. (|TED ) . Contributions from DIS from the pionic fields themselves, which are needed to 



balance overall momentum conservation, were considered in Refs. p8| , |30| - p^ , p^1 , however, 
their role is most evident only at small x (x < 0.2). 

The bulk of the suppression of the EMC ratio (0) at x ~ 0.6 in the binding model 
can be attributed to the fact that the average value for the interacting nucleon light-cone 
fraction is less than unity. For the case of the deuteron, this corresponds to the average 
spectator light-cone fraction (a s ) > 1, which is contrary to what one would have from Fermi 
motion alone, where the average a s is < 1. A relatively minor role is played by the structure 
function of the bound nucleon itself — the only requirement is that it be a monotonically 
decreasing function of x |]8|, |28|j35[| . This is clearly the case for the on-shell structure function, 



and since the off-shell behavior of the bound nucleon structure function is unknown, most 
early versions of the binding model simply neglected the possible dependence on p 2 with 
the expectation that it is not large for weakly bound systems. Only very recently has the 
issue of off-shell dependence in the bound nucleon structure function been addressed P1,|T0 



where the first attempts to construct models for the p 2 dependence of F^lq were made. 
Note that a consequence of assuming the absence of any off-shell effects in F^J is that the 
tagged structure function ratio G in Eq. (|5B|) , normalized to the corresponding ratio for a 
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free proton, would be unity (see Fig.5 below). Any observed deviation of this ratio from 
unity would therefore be a signal of the presence of nucleon off-shell effects. 

If in a dilute system such as the deuteron the nucleon off-shell effects do not play a major 
role (at least at x ^ 0.7), one could expand the effective nucleon structure function in a 
Taylor series about p 2 = M 2 : 



F eff \ 



x,p 2 ,Q 2 ) = F 2N (x,Q 2 ) + (p 2 -M 2 



dF^ f (x,p 2 ,Q 2 ) 



dp 2 



+ 



(34) 



=M 2 



Here the off-shell dependence is determined, to order p 2 /M 2 , entirely by the derivative 
of F 2 J with respect to p 2 . In order to model this correction a microscopic model of the 



nucleon structure is required P,ID,r^|. In any generic quark-parton model, the effective 
nucleon structure function can be written as an integral over the quark momentum p q of the 
quark spectral function p: 



( x , P J Q 



j d P 2 g p(p 2 q,P 2 ,X,Q 2 



(35) 



To proceed from Eq.(|35|) requires additional assumptions about the quark spectral function. 
The simplest is to assume that the p 2 and p 2 dependence in p is factored [H],^], which then 

leads to an explicit constraint on the p 2 derivative of from baryon number conservation 
in the deuteron: 



i dx dF^ f (x,p 2 ,Q 2 



X 



dp 2 



(36) 



=M 2 



[x,p 2 ,Q 2 )/dp 2 to be determined from the x-dependence of the on- 



This then allows dF 2 f J ( 
shell structure function F 2 n(x,Q 2 ) in terms of a single parameter, the squared mass of the 
intermediate state "diquark" system that is spectator to the deep-inelastic collision, (p—p q ) 2 - 
One can obtain a good fit to the on-shell nucleon structure function data in terms of this 
model if one restricts the spectator "diquark" mass to be in the range (p — p q ) 2 ~ 2 — 4 
GeV 2 m. 



A more microscopic model which does not rely on the factorization of the p 2 and p 2 de- 
pendence in p was discussed in Ref. [[|. The quark spectral function there was determined 
entirely from the dynamics contained in the nucleon-quark-spectator "diquark" vertex func- 
tion, r(p,p q ). Within the approximation discussed in Section [1TA] , taking the positive energy 
nucleon projection only, 



P(pI,P 2 ,X,Q 2 



Tr 



+ M) T{p,p q ) 



9 q - m q 



,-i 



(37) 



where m q is the quark mass. Angular momentum conservation allows two forms for the 



vertex function T, namely scalar and pseudo-vector. In Refs. 0,|Tl],|l4[] it was found that 
the on-shell data could be well described in terms of only a few of the many possible Dirac 
structures for V. In particular, the vertices were chosen to be cx / and 7^75. The momentum 
dependence of the vertex functions, on the other hand, is more difficult to derive, and must 
in practice be either parameterized or calculated by solving bound state Faddeev equations 
in simple models of the nucleon ||36|| . In order to obtain realistic static properties of the 
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nucleoli, and to account for the bound nature of the nucleon state, the vertex functions 
must have the form: 



r(p,p 9 ) oc 



(™> 2 q ~ Pg) 

(a 2 - p, 



(38) 



where the parameters A and n are fixed by comparing with the quark distribution data, and 
the overall normalization is fixed by the baryon number conservation condition for both the 
nucleon and deuteron structure functions In the comparisons in Sec. HID we use 

the parameters from the analysis of Ref . [H] . 



B. Color Screening Model of Suppression of Point-like Configurations in Bound 

Nucleons 



A significant EMC effect in inclusive (e, e') reactions occurs for x ~ 0.5-0.6 which corre- 
sponds to the high-momentum component of the quark distribution in the nucleon. There- 
fore the EMC effect in this x range is sensitive to a rather rare component of the nucleon 



wave function where 3 quarks are likely to be close together [17,25]. It is assumed in this 
model that for large x the dominant contribution to F2n{x,Q 2 ) is given by the point-like 
configurations (PLC) of partons which weakly interact with the other nucleons. Note that 
due to scaling violation F 2 n(x,Q 2 ) at x ^ 0.6, Q 2 ^ 10 GeV 2 , is determined by the non- 
perturbative nucleon wave function at x > 0.7. Thus it is actually assumed that in the 
nonperturbative nucleon wave function point-like configurations dominate at x ^ 0.7. The 
suppression of this component in a bound nucleon is assumed to be the main source of the 
EMC effect in inclusive deep-inelastic scattering P7[ , |25| . Note that this suppression does 
not lead to a noticeable change in the average characteristics of nucleons in nuclei |25] . 

To calculate the change of the probability of a PLC in a bound nucleon, one can use a 
perturbation series over a small parameter, k, which controls corrections to the description 
of a nucleus as a system of undeformed nucleons. This parameter is taken to be the ratio of 
the characteristic energies for nucleons and nuclei: 



k=| (U A )/AE A |~1/10, 



(39) 



where (U A ) is the average potential energy per nucleon, (U A ) U»i~ — 40 MeV, and AE A ~ 
M* — M ~ 0.5 GeV is the typical energy for nucleon excitations within the nucleus. Note 
that AEd ~ 2(Ma — M) ~ 0.6 GeV, since the N — A component in the deuteron wave 
function is forbidden due to the zero isospin of the deuteron. 

To estimate the deformation of the bound nucleon wave function we consider a model 
where the interaction between nucleons is described by a Schrodinger equation with potential 
V(Rij,yi,yj) which depends both on the inter-nucleon distances (spin and isospin of nucle- 
ons) and the inner variables yi and yj, where y% characterizes the quark-gluon configuration 
in the i-th nucleon [ 17 ,25|,37]. The Schrodinger equation can be represented as: 



2m 



- E v 2 + E viRiv + E H o(vi) 

N 4 i A 4 



i,0 



rj}{y u Rij) = Etp(y i ,R ij ). 



(40) 
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Here Ho(yi) is the Hamiltonian of a free nucleoli. In the nonrelativistic theory of the nucleus 
the inter-nucleon interaction is averaged over all yi. Thus the nonrelativistic U(Rij) is related 
to V as: 

u ( R ij)= Yl (<M^)<M%-) \V(R ij ,y h y j: y h y j ) | 0jv(|/i)0iv (%■)}, ( 41 ) 
yi,yj,yi,yj 

where 4>N{Vi) is the free nucleon wave function. The unperturbated wave function is the 
solution of Eq. (|40D with the potential V replaced by U. Treating {U — V)/(Ei — En) as a 
small parameter, where Ei is the energy of an intermediate excited nucleon state, one can 
calculate the dependence of the probability to find a nucleon in a PLC to the momentum of 
a nucleon inside the nucleus. One finds that this probability is suppressed as compared to 
the similar probability for a free nucleon by the factor f2T 



S A (k 2 ) « 1 - 4(k 2 /2M + e a )/AE A , (42) 

where AE A = (Ei — En) ~ M* — M, in first order of the perturbation series. An estimate 
of higher order terms gives AFT 



Mk 2 ) = (1 + z)-\ z = (k 2 /M + 2e A )/AE A . (43) 

The x dependence of the suppression effect is based on the assumption that the PLC con- 
tribution in the nucleon wave function is negligible at x < 0.3, and gives the dominant 
contribution at x > 0.5 p5|,|38| . We use a simple linear fit to describe the x dependence 



between these two values of x f3~8fl . One can then obtain an estimate for R A in Eq.(|2|) for 
large A at x ~ 0.5, 

Ra(x) Uo.5~ S A (k 2 ) « 1 + ~ 0.7 - 0.8, (44) 

since here Fermi motion effects are small. The excitation energy AE A for the compressed 
configuration is estimated as AE A ~ (M(1400) — M) — (M(1680) — M) ~ (0.5-0.8) GeV, 
while (U A ) ~ —40 MeV. Since (U A ) ~ (pa(j)) the model predicts also the A dependence of 
the EMC effect, which is consistent with the data [ |T7|j . 

For the deuteron target we can deduce from Eq . (f44l) using (Ud) / (UFe) ~ 1/5, 

Rd(x,Q 2 ) 1^0.5- 0.94 - 0.96. (45) 

This number may be somewhat overestimated because, as discussed above, due to the 
isoscalarity of D, low-energy excitations in the two-nucleon system are forbidden, leading to 
AE D > AE A . 

This model represents one of the extreme possibilities that the EMC effect is solely the 
result of deformation of the wave function of bound nucleons, without attributing any extra 
momentum to be carried by mesons. A distinctive feature of this explanation of the EMC 
effect is that the deformation of the nucleon should vary with inter-nucleon distances in 
nuclei (with nucleon momentum in the nucleus). For |k| ~ 0.3 — 0.4 GeV/c the deviation 
from the conventional quantum-mechanical model of a deuteron is expected to be quite large 
(factor ~ 2). Actually, the size of the effects may depend not on k 2 only but on and 
a s separately, because the deformation of a bound nucleon may be more complicated than 
suggested by this simple model. 
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C. QCD Radiation, Quark Derealization. 



It was observed in Refs. p9|^0| that the original EMC data could be roughly fitted as: 

jF 2A (x, Q 2 ) = l -F 2D (x, Q 2 UQ 2 )) , (46) 

with £,Fe{Q 2 ) ~ 2 for Q 2 w 20 GeV 2 , the so-called dynamical rescaling. The phenomeno- 
logical observation has been interpreted as an indication that gluon radiation occurs more 
efficiently in a nucleus than in a free nucleon (at the same Q 2 ) due to quark derealization, 
either in a bound nucleon (or in two nearby nucleons) [09| ,p|-^3| or in the nucleus as a whole 



H 



The Q 2 dependence of £(Q 2 ) follows from the requirement that both sides of Eq . (|46|) 
should satisfy the QCD evolution equations. In the leading logarithmic approximation one 
has: 

UQ 2 ) = u{QlT QCD(Ql)/aQCDm - (47) 

Experimentally dlnF2D{x, Q 2 )/dlnQ 2 is positive if x > xq and negative if x < xq, where 
xq = 0.15 ± 0.05. So Eq.(f4~6"D predicts that the EMC effect should vanish at x ~ xq. 

If the confinement size in a nucleus, A^, is larger than that in a free nucleon, one may 
expect that the Q 2 evolution of the parton distributions in nuclei (the bremsstrahlung of 
gluons and quarks) may start at Ql(A) < Ql(N). To reproduce Eq.(|46|) one should have: 

Ql(A)/Ql(N) = ml(A))}- 1 . (48) 

Assuming on dimensional grounds that the radii of quark localization in a nucleus, A^, and 
in a nucleon, X N , are related via: 

Ql(A)\\ = Ql(N)X 2 N (49) 
to the scale for the onset of evolution, Qq, one finally obtains: 

\ A /\ N ^aQUA)) 1/2 . (50) 

A fit to the original EMC data using Eq.(|§ leads to £ Fe (20 GeV 2 ) =2. For fx 2 Fe = 0.67 GeV 2 
and Ajjg rs 250 MeV this corresponds to: 

\ Fe /\ D « 1.15 A^ 200 /A D w 1.19. (51) 

Since in this model the delocalization is approximately proportional to the nuclear den- 
sity, one can expect that the effect is also proportional to the k 2 of a bound nucleon. Fixing 
the parameters of the model to fit the Fe data, and assuming X(k)/X = 1 + ak 2 we obtain 
a w 0.4/ < k 2 >Fe, where < k 2 >Fe~ 0.08 GeV 2 /c 2 . Using this expression for X(k) one 
can calculate the id{Q 2 ,k) according to Eqs. (f£7|)-(|50|) and express the effective structure 
function of a bound nucleon as: 

F^ f (x,a s ,p T ,Q 2 ) = F 2N (x,Q 2 £(Q 2 ,k)) , (52) 

where k is defined through a s and according to Eq. (|23|) . 
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D. Numerical estimates 



Comparison of predictions of the models for the nuclear EMC effect considered in the 
preceding Sections is most meaningful in the kinematic range where, firstly, FSI effects are 
small, and secondly, the instant-form and light-cone prescriptions for the deuteron spectral 
function lead to similar results. 

Direct calculation of the FSI contribution to the cross section would require knowledge of 
the full dynamics of the final N-X system, which is a practically impossible task given the 
present level of understanding of nonperturbative QCD. However, it is possible to estimate 
the uncertainty which would be introduced through neglect of FSI, by using the calculation 
of FSI effects in the high-energy d(e, e'p)n (break-up) reaction in Refs. |23| , |24|| , and replacing 
the p-n rescattering cross section by an effective cross section for the p(n)-X interaction, 
Eq. (p7|) . For the effective cross section a e // one can use the results of the recent analysis |21| 



of soft neutron production in the high-energy deep-inelastic scattering of muons from heavy 
nuclei J22|, which yielded an upper limit of a e ff ~ 20 mb. To be on the conservative side, 
in the following estimates we therefore use the value of a e ff ~ 20 mb suggested by string 
models of FSI (for a recent summary see Ref. [fHJ). Furthermore, by retaining only the 
imaginary part of the spectator nucleon rescattering amplitude, one obtains an upper limit 
of the FSI effect, since the real part will contribute to elastic rescattering only, effectively 
suppressing the value of cr e ff- In Fig.l we illustrate the a s dependence of the ratio of 
the (light-cone) spectral function including FSI effects within the DWIA, Eq.(^Tj), to that 
calculated without FSI effects. At extreme backward kinematics (pt ~ 0) one sees that 
FSI effects contribute less than ~ 5% to the overall uncertainty of the d(e, e'N)X cross 
section for a s < 1.5. As mentioned above, this number can be considered rather as an upper 
limit on the uncertainties due to FSI. At larger p T (> 0.3 GeV/c), and small a s (~ 1), 
the double scattering contribution (which is not present for the extreme backward case in 
Eq.([27D) plays a more important role in FSI [p3fl . Because its sign is positive, it tends to 
cancel some of the absorption effects of FSI at large pt (for a detailed discussion of the 
double scattering contribution in FSI see Ref. [p3| ). Note also that the FSI effects do not 
change significantly with energy for fixed Q 2 . Thus, for the ratios discussed, where the cross 
sections are compared for the same Q 2 but different x, the changes due to FSI effects are 
even smaller. 

To extract unambiguous information from the semi-inclusive cross section ratios on the 
medium modifications of the nucleon structure discussed in this Section requires one to 
establish the regions of kinematics where the differences between the various prescriptions for 
the deuteron spectral function are minimal. In Fig. 2 we illustrate the a s and px dependence 
of the ratio of spectral functions calculated in the light-cone ( |2"2"D and instant-form ( |TTD 
approaches. For pt < 0.1 GeV/c the light-cone and instant-form predictions differ up to the 
20% for the entire range of a s < 1.5. However, choosing the isolated values of a s < 1.2 or 
a s ~ 1.4 one can confine the uncertainty in the spectral function to within 10%, which will 
offer the optimal conditions in which to study the nucleon structure modification. 

In Fig. 3 we show the ot s dependence of the ratio of the effective proton structure function, 
, for extreme backward kinematics, pr = 0, to the structure function of a free proton. 
As expected, the suppression of the ratio in the version of the binding model with explicit 
nucleon off-shell corrections is quite small, ^ 10% for a s < 1.5, reflecting the relatively minor 
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role played here by off-shell effects in the nucleon structure function. (Note that in versions 
of the binding model in which there is no p 2 dependence in the effective nucleon structure 
function this ratio would be unity.) The effects in the PLC suppression and rescaling models 
in Fig. 3 are somewhat larger. For a neutron target one predicts similar results, however, 
the neutron structure function is not as well known experimentally as the proton due to the 
absence of free neutron targets (see Section IV). For this reason we restrict our discussion 
to ratios of proton structure functions only. 

In Fig. 4 we illustrate the dependence of F 2 l /F%p on the variable (p 2 — 2Me)/M 2 , where 
p 2 is the bound nucleon virtuality defined in Eq.(pH|) and e = —2.2 MeV is the deuteron 
binding energy. The comparison in Fig. 4 is done for different values of x and a s . It is 
noticeable that because of the negligible amount of the PLC component in the nucleon wave 
function at x ^ 0.3, the PLC suppression model predicts no modification of the structure 
function in this region. On the other hand, at x = 0.6 it predicts a maximal effect because 
of the PLC dominance in the nucleon wave function here. 

Because the momentum (virtuality) dependent density effect generates the modification 
of the bound nucleon structure functions in the PLC suppression and rescaling models, the 
ratios for these models at x > 0.5-0.6 vary similarly with a s and (p 2 — 2Me)/M 2 , as does 
also the off-shell model ratio. However the mechanism for such a variation is different, which 
is clearly seen in Fig. 5, where the x dependence of the same ratio is represented at different 
values of a and fixed pt = 0. The curves for the off-shell model extend only to a; ~ 0.7 
because for larger x values the approximations discussed in Sec . [II A 1| involved in obtaining 
Eq. fllBD become numerically less justified [P,PH|. 

To further reduce any uncertainties due to the deuteron spectral function in the model 
comparisons, we concentrate on the predictions for the ratio G(a s ,pT, Xi, x 2 , Q 2 ), defined in 
Eq. (|3"3"D , of experimentally measured cross sections at two different values of x |25| . Since the 
function G is defined by the ratio of cross sections at the same a s and pr, any uncertainties 
in the spectral function cancel. This allows one to extend this ratio to larger values of a s , 
thereby increasing the utility of the semi-inclusive reactions when analyzed in terms of this 
function. Figure 6 shows the a s and Q 2 dependence of G(a s ,pT, x\, x 2 , Q 2 ) at p T = 0. The 
values of X\ and x 2 are selected to fulfill the condition £i/(2 — a s ) = 0.6 (large EMC effect 
in inclusive measurements) and x 2 /(2 — a s ) = 0.2 (essentially no EMC effect in inclusive 
measurements). Again, the PLC suppression and Q 2 rescaling models predict a much faster 
drop with a s than does the binding/off-shell model, where the a s dependence is quite weak. 



IV. EXTRACTION OF THE NEUTRON/PROTON STRUCTURE FUNCTION 

RATIO 

The presence of an EMC effect in the deuteron leads to substantial suppression of the 
deuteron structure function at large x compared to what one would expect from models 
without non-nucleonic degrees of freedom. For example, the estimate of Ref. ||25|| , 

F 2D (x,Q 2 ) 1 2F 2A (x,Q 2 



F 2p (x, Q 2 ) + F 2n (x, Q 2 ) 4 AF 2D (x, Q 2 



(53) 

A~60;0.3<x<0.7 



which is valid for a rather wide class of models in which the EMC effect is proportional to 
the mean value of p 2 in nuclei, gives a ratio ~ 3-5% below unity at x ~ 0.6-0.7. Calculations 
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of the deuteron structure function in models in which binding effects are explicitly taken 



into account also produce similar effects [1 1 



Inclusion of the EMC effect in the extraction of the the neutron structure function from 
the inclusive eD scattering data leads to significantly larger values for F 2n /F 2p than the 1/4 
value obtained in analyses in which only Fermi motion is included. The values for F 2n jF 2p 
with inclusion of the EMC effect at x ~ 0.6 f45"l-|47] are in fact much closer to the expectation 
of 3/7 from perturbative QCD, predicted by Farrar and Jackson [58]. Therefore observation 
of a value of F 2n /F 2p higher than the 1/4 extracted from inclusive data in the early analyses 
would by itself serve as another proof of the presence of the EMC effect in the deuteron. 

Although other methods to obtain the large- x n/p ratio (or the d/u ratio) have been 
suggested, none has so far been able to clearly discriminate between the different x — > 1 
limits for F 2n /F 2p (namely, d/u — > 0, which is the minimal possible value allowed in the 
parton model, which corresponds to F 2n /F 2p — > 1/4, and d/u — > 1/5 in perturbative QCD 
4"8f). For example, with v and V beams on proton targets one can in principle measure the u 



and d quark distributions separately, however, the statistics in v experiments in general are 
relatively poor. Another possibility would be to extract the d/u ratio from charged lepton 



asymmetries at large rapidities, in P^-boson production in pp scattering |49| , although here 
it may be some time still before a sufficient quantity of large-rapidity events at the CDF at 
Fermilab are accumulated. 

On the other hand, with tagged deuteron experiments planned for HERMES, a study 
of the tagged structure functions may allow a resolution of this ambiguity H§,^(J. It is 
important that the ratio of tagged structure functions interpolated to the nucleon pole 
should be exactly equal to the free nucleon ratio — this is the analog of the Chew-Low 
interpolation for the pion case: 



F 2n (x,Q 2 ) _ F 2 e l f (x/(2-a s ),p T ,Q 2 
F 2p (x,Q 2 ) ' 



F^ f (x/(2-a s ),p T ,Q 2 ) 



(54) 



a s !»l,PT~0 



In practice the data cannot be accumulated for too small p. However, we observed above 
that deviations of the ratio from the free limit is proportional to p 2 with a good accuracy. 
Hence, if one samples the data as a function of p 2 interpolation to the pole p 2 — M 2 = 2Me 
should be smooth. In practice, considering momentum intervals of 100-200 MeV/c and 
200-350 MeV/c would be sufficient. A potential problem with Eq . (|54]) is that at very large 
x (x > 0.7) the factorization approximation itself breaks down and higher order corrections 
to Eq. ([IB]) , which are oc p 4 , must be included if one wants accuracy to within a few %. To 
be on the safe side one should therefore restrict the analysis to smaller spectator momenta, 
below 200 MeV/c. 



V. CONCLUSION 

Despite the many years of study of the deviations from unity of the ratios of nuclear to 
deuteron cross sections in inclusive high-energy scattering, we have been unable to isolate the 
microscopic origin of the nucleon structure modification in the nuclear environment |IJ-§[ - 
the effect can be described in terms of a number of models based on quite disparate physical 
assumptions. In this paper we have argued that by allowing greater accessibility to kinematic 
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variables not available in inclusive reactions, semi-inclusive deep-inelastic processes offer the 
possibility to make further progress in understanding the origin of the nuclear EMC effect. 

In particular, measurements of tagged structure functions of the deuteron can probe the 
extent of the deformation of the intrinsic structure of a bound nucleon. Taking ratios of 
semi-inclusive cross sections at different values of x enables the cancellation in the tagged 
structure function ratio, Eq . (|33|) , of the dependence on the deuteron wave function, thus 
permitting the nucleon structure to be probed directly. Our results show that possible 
contamination of the signal due to the final state interactions of the spectator nucleon with 
the hadronic debris can be minimized by tagging only on the slow backward nucleons in the 
target fragmentation region. This may then allow one to discriminate between models in 
which the EMC effect is attributed entirely to the nucleon wave function deformation, and 
ones in which the effect arises from more traditional descriptions in terms of meson-nucleon 
degrees of freedom associated with nuclear binding. 

As a by-product of the semi-inclusive measurements, one may also be able to extract 
information on the large-x behavior of the neutron to proton structure function ratio, by 
detecting recoiling protons and neutrons with small transverse momentum in the extreme 
backward kinematics. Extracting this ratio from inclusive ep and eD data is fraught with 
large uncertainties arising from different treatments of the nuclear physics in the deuteron. 
Observation of an asymptotic value for F 2n /F2 P which is larger than the 'canonical' 1/4 
would in itself be proof of the presence of an EMC effect in the deuteron. 

First information about the spectra of backward nucleons in eD scattering is likely to 
come from the Jefferson Lab experiment # 94-102 0, and from the HERMES experiment at 
HERA, where the necessary counting rate may be achieved after an upgrade of the detector 



5j| . Having two energy ranges would be very useful for checking the basic production 
mechanism, and understanding backgrounds and corrections due to final state interactions, 
which are likely to depend substantially on the incident energy. Final state interaction effects 
can also be tracked by studying the production of nucleons from heavier nuclei. Probing the 
quark-gluon structure of short-range correlations with heavy nuclei targets should further 
enable one to determine whether nucleon deformations predominantly depend on the nucleon 
momentum, or also on A. 
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FIGURES 



FIG. 1. The a s dependence of the ratio of cross sections calculated with FSI effects within the 
DWIA, and without FSI effects. The curves correspond to different values of the spectator nucleon 
transverse momentum (in GeV/c). 

FIG. 2. The a s dependence of the ratio of cross sections calculated within the light-cone and 
instant-form approaches. The curves correspond to different values of the spectator nucleon trans- 
verse momentum (in GeV/c). 

FIG. 3. The a s dependence of /F 2p for x = 0.6 and p T = 0. Dashed line is the PLC 
suppression model, dotted is the rescaling model, and dot-dashed the binding/off-shell model. 

FIG. 4. The (p 2 - 2Me)/M 2 dependence of F^ f /F 2p , for a s = 1.2 and 1.4, and x = 0.3 and 
0.6. Curves are as in Fig.3. 

FIG. 5. The x dependence of F^ p / ' F 2p for a s = 1.2 and 1.4, with px = 0. Curves are as in Fig.3. 

FIG. 6. The a s dependence of G(a s ,pT, x\, x 2 , Q 2 ), with x\ = x/{2 — a s ) = 0.6 and 
X2 = x/(2 — a s ) = 0.2, for pt = 0. G e ^ (a s ,pT,xi, x 2 ,Q 2 ) is normalized to G e ^ (a s ,pT, x\,x 2 , Q 2 ) 
calculated with the free nucleon structure function. Curves are as in Fig.3. 
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